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We calculate the free energy of a hot gas of electrons and photons to three loops using the hard- 
thermal-loop perturbation theory reorganization of finite-temperature perturbation theory. We 
calculate the free energy through three loops by expanding in a power series in mo/T, rrif/T, and 
e 2 , where mo and m/ are thermal masses and e is the coupling constant. We demonstrate that the 
hard-thermal-loop perturbation reorganization improves the convergence of the successive approxi- 
mations to the QED free energy at large coupling, e ~ 2. The reorganization is gauge invariant by 
construction, and due to cancellation among various contributions, we obtain a completely analytic 
result for the resummed thermodynamic potential at three loops. Finally, we compare our result 
with similar calculations that use the ^-derivable approach. 
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I. INTRODUCTION 

The thermodynamic functions for hot field theories can 
be calculated as a power series in the coupling constant g 
at weak coupling. One is primarily interested in calculat- 
ing the free energy from which the pressure, energy den- 
sity, and entropy can be obtained using standard thermo- 
dynamic relations. In the early 1990s the free energy was 
calculated to order o 4 in Refs. [l], Q for massless scalar 
4> 4 theory, in Ref. p[ for QED and in Ref. Q for non- 
Abclian gauge theories. The corresponding calculations 
to order g 5 were carried out in Refs. [J, Refs. @, 0] 
and Refs. [1, Q, respectively. Recent results have ex- 
tended the calculation of the QCD free energy by deter- 
mining the coefficient of the g 6 log(g) contribution [l(J. 
For massless scalar the perturbative free energy is now 
known to order g 6 [111 an d <7 8 l°g(ff) 03 ■ 

Unfortunately, the resulting weak-coupling approxima- 
tions, truncated order-by-order in the coupling constant, 
are poorly convergent unless the coupling constant is ex- 
tremely small. For example, simply comparing the mag- 
nitude of low-order contributions to the Nf = 3 QCD 
free energy one finds that the g^ contribution is smaller 
than the g^ contribution only for g s < 0.9 (a s < 0.07). 
This is a troubling situation since at phenomenologically 
accessible temperatures near the critical temperature for 
the QCD deconfinement phase transition, the strong cou- 
pling constant is on the order of g s ~ 2. We therefore 
need methods which can provide reliable approximations 
to QCD thermodynamics at intermediate coupling. 

The poor convergence of finite-temperature perturba- 



tive expansions of the free energy is not limited to QCD. 
The same behavior can be seen in weak- coupling ex- 
pansions in scalar field theory [l| and QED [y, @| . In 
Fig. [1] we show the successive perturbative approxima- 
tions to the QED free energy. As can be seen from this 
figure, at couplings larger than e ~ 1 the QED weak- 
coupling approximations also exhibit poor convergence. 
For this reason a concerted effort has been put forth to 
find a reorganization of finite-temperature perturbation 
theory which converges at phenomenologically relevant 
couplings. Here we will focus on the QED free energy 

There arc several ways of systematically reorganiz- 
ing the perturbative expansion to improve its conver- 
gence and the various approaches have been reviewed 
in Refs. [3, EH, [IB]- Here we will focus on the 
hard-thermal- loop perturbation theory (HTLpt) method 
[13, [H, [H, [M EJ. The HTLpt method is inspired by 
variational perturbation theory [H, H, H, SH] . HTLpt 
is a gauge-invariant extension of screened perturbation 
theory (SPT) Qj|, [2y, [27|, [28| , which is a perturbative re- 
organization for finite-temperature massless scalar field 
theory. In the SPT approach, one introduces a single 
variational parameter which has a simple interpretation 
as a thermal mass. In SPT a mass term is added to and 
subtracted from the scalar Lagrangian, with the added 
piece kept as part of the free Lagrangian and the sub- 
tracted piece associated with the interactions. The mass 
parameter is then required to satisfy a variational equa- 
tion which is obtained by the principle of minimal sensi- 
tivity. 

This naturally led to the idea that one could apply a 
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FIG. 1: Successive perturbative approximations to the QED 
pressure (negative of the free energy). Each band corresponds 
to a truncated weak-coupling expansion accurate to order e 2 , 
e 3 , e 4 , and e 5 , respectively. Shaded bands correspond to vari- 
ation of the renormalization scale /j, between tvT and 4ttT. 



similar technique to gauge theories by adding and sub- 
tracting a mass in the Lagrangian. However, in gauge 
theories, one cannot simply add and subtract a local 
mass term since this would violate gauge invariance. In- 
stead one adds and subtracts to the Lagrangian a hard- 
thermal-loop (HTL) improvement term. The free part of 
the Lagrangian then includes the HTL self-energies and 
the remaining terms are treated as perturbations. Hard- 
thermal-loop perturbation theory is a manifestly gauge- 
invariant approach that can be applied to static as well 
as dynamic quantities. SPT and HTLpt have been ap- 
plied to four and two loops [H, 03, Gj, Ql, EI EH E| , 
respectively, and convergence is improved compared to 
the weak-coupling expansion. 

In this paper we calculate the pressure in QED to 
three-loop order in HTLpt. This will set the stage for 
the corresponding calculation in full QCD. We deter- 
mine leading-order (LO), next-to-leading-order (NLO), 
and next-to-next-to-leading-ordcr (NNLO) expressions 
for the HTLpt pressure. At NNLO the expression is 
entirely analytic and gives a well-defined gap equation 
(variational equation) for the electron and photon screen- 
ing masses. As we will show, the NLO and NNLO 
HTLpt rcsummcd QED free energy give approximations 
which show improved convergence for couplings as large 
as e ~ 2.5 (see Fig. [7]). In addition, we compare our 
results to those obtained using the 2PI ^-derivable ap- 
proach [3l|, EI| and show that at three loops the agree- 
ment between the HTLpt and <I>-derivable approaches is 
quite good. 

The structure of the paper is as follows. We give a 
brief summary of HTLpt in Sec. II. In Sec. Ill, we list the 
expressions for the one-, two-, and three- loop diagrams 
that contribute to the thermodynamic potential. In Sec. 



IV, we expand the sum-integrals in the mass parameters, 
and in Sec. V, the free energy is calculated. We summa- 
rize and draw our conclusions in Sec. VI. In Appendix 
A, we give the Feynman rules for HTLpt in Minkowski 
space. In Appendixes B and C, we list all sum-integrals 
and integrals needed in the calculations. 



II. HTL PERTURBATION THEORY 

The Lagrangian density for massless QED in 
Minkowski space is 



£qed = -\f„ 



£ g f + £ g h + A£q ED 



(1) 



Here the field strength is = d»A v - d v A» and the 
covariant derivative is = + ieA^. The ghost term 
£ g h depends on the gauge- fixing term £ g f . In this paper 
we choose the class of covariant gauges where the gauge- 
fixing term is 



(2) 



with £ being the gauge-fixing parameter. In this class of 
gauges, the ghost term decouples from the other fields. 

The perturbative expansion in powers of e generates 
ultraviolet divergences. The renormalizability of pertur- 
bative QED guarantees that all divergences in physical 
quantities can be removed by renormalization of the cou- 
pling constant a = e 2 /4ir. There is no need for wave- 
function renormalization, because physical quantities are 
independent of the normalization of the field. There is 
also no need for renormalization of the gauge parameter, 
because physical quantities are independent of the gauge 
parameter. 

Hard-thcrmal-loop perturbation theory is a reorganiza- 
tion of the perturbation series for thermal gauge theories. 
In the case of QED, the Lagrangian density is written as 



£ — (£qed + £htl) 



AC 



HTL 



(3) 



The HTL improvement term is 



£htl 



--(l-<5)my^ a 



(yd) 2 / y 



+ (1 - S) imjip^ 



y-D 



0, 



(4) 



The 



where y M = (l,y) is a light-like four- vector, and 
represents an average over the directions of y. 
term (Q| has the form of the effective Lagrangian that 
would be induced by a rotationally-invariant ensemble 
of charged sources with infinitely high momentum. The 
parameter mjj can be identified with the Debye screen- 
ing mass and the parameter m,f can be identified as the 



3 



induced finite-temperature electron mass. HTLpt is de- 
fined by treating 5 as a formal expansion parameter. 

The HTL perturbation expansion generates ultravio- 
let divergences. In QED perturbation theory, renormal- 
izability constrains the ultraviolet divergences to have 
a form that can be cancelled by the counterterm La- 
grangian A£qed- We will demonstrate that renormal- 
ized perturbation theory can be implemented by includ- 
ing a counterterm Lagrangian A£htl among the inter- 
action terms in There is no proof that the HTL 
perturbation expansion is renormalizable, so the general 
structure of the ultraviolet divergences is not known; 
however, it was shown in previous papers [20l . [2l| that 
it was possible to renormalize the next-to-leading-order 
HTLpt prediction for the free energy of QED using only 
a vacuum counterterm, a Debye mass counterterm, and 
a fcrmion mass counterterm. In this paper we will show 
that renormalization is also possible at NNLO. 

The counterterms necessary are 
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(5) 
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Physical observables are calculated in HTLpt by ex- 
panding them in powers of 5, truncating at some specified 
order, and then setting 5=1. This defines a reorgani- 
zation of the perturbation series in which the effects of 
the m 2 D and m 2 terms in ((4]) are included to all orders 
but then systematically subtracted out at higher orders 
in perturbation theory by the 5m 2 D and 5m 2 terms in Q . 
If we set 8=1, the Lagrangian ([3]) reduces to the QED 
Lagrangian |T]). 

If the expansion in 8 could be calculated to all or- 
ders, the final result would not depend on tod or to/ 
when we set 8 = 1. However, any truncation of the ex- 
pansion in 6 produces results that depend on mo and 
to/. Some prescription is required to determine tod and 
to/ as a function of T and a. We choose to treat both 
as variational parameters that should be determined by 
minimizing the free energy. If we denote the free energy 
truncated at some order in 5 by f2(T, a, tod, to/, /i, 5), 
our prescription is 



d 



d 



dm 



f 



fi(T, a, tod, to/, /x, 5 = 1) = 0, (9) 
fi(T, a, tod, to/, /x, 5 = 1) = 0. (10) 



Since £l(T, a, tod, To/, fx, 5 = 1) is a function of the vari- 
ational parameters tod and to,/ , we will refer to it as the 
thermodynamic potential. We will refer to the variational 



equations © and (fTO)) as the gap equations. The free en- 
ergy T is obtained by evaluating the thermodynamic po- 
tential at the solution to the gap equations ((9]) and (fTO]) . 
Other thermodynamic functions can then be obtained by 
taking appropriate derivatives of T with respect to T. 



III. DIAGRAMS FOR THE THERMODYNAMIC 
POTENTIAL 



In this section, we list the expressions for the diagrams 
that contribute to the thermodynamic potential through 
order 5 2 in HTL perturbation theory. The diagrams are 
shown in Figs. [21 [3] and [|] Because of our dual truncation 
in tod, to/, and e the diagrams listed in Fig. 3] do not 
contribute to our final expression so we will not explicitly 
list their integral representations. The expressions here 
will be given in Euclidean space; however, in Appendix 
lAl we present the HTLpt Feynman rules in Minkowski 
space. 

The thermodynamic potential at leading order in HTL 
perturbation theory for QED with Nf massless electrons 



^lo = ^la + N f T xh + A £ • (11) 
Here, T\ a is the contribution from the photons 

Fi* = -\f p ^ d - *) lo § 1- A t(P)} + log A L (P)} . (12) 

The transverse and longitudinal HTL propagators At(P) 
and Al(P) are given in (jMTj) and (|A52j) . The electron 
contribution is 



Fib = -i logdet[f-E(P)] 
T{p} 



(13) 



The leading-order vacuum counterterm AqSq is given by 

1 



A £ = 
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(14) 



The thermodynamic potential at next-to-leading order 
(NLO) in HTL perturbation theory can be written as 



^NLO = ttLO+A r /(^2a+^2b+^2d)+^2c + Al£ 

d 
dm 



d d 
+A 1 m 2 D —^-n LO + Aim 2 — jJ^lq , (15) 

0Tft j-y OTUj 



where A\£q, Aim 2 D , and A\m 2 are the terms of order 5 
in the vacuum energy density and mass counterterms: 
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FIG. 2: Diagrams contributing through NLO in HTLpt. The undulating lines are photon propagators and the solid lines are 
fermion propagators. A circle with a II indicates a photon self-energy insertion and a circle with a E indicates a fermion 
self-energy insertion. All propagators and vertices shown are HTL-resummed propagators and vertices. 
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FIG. 3: Diagrams contributing to NNLO in HTLpt which contribute through order e 5 . The undulating lines are photon 
propagators and the solid lines are fermion propagators. A circle with a II indicates a photon self-energy insertion and a circle 
with a E indicates a fermion self-energy insertion. The propagators are HTL-resummed propagators and the black dots indicate 
HTL-resummed vertices. The lettered vertices indicate that only the HTL correction is included. 



The contributions from the two-loop diagrams with 
electron-photon three- and four-point vertices are 



•^2a = \e 



2 JP{Q} 



Tr rr"(P, Q, R)S(Q)T V (P, Q, R)S(R)] 



xA""(?) , 
[T^(P,-P. 

xA^(P) , 



•^2b = \e 2 4 Tr[T^(P,-P,Q,Q)S(Q)] 
2 JP{Q} 



(19) 



(20) 



where R = Q — P. 

The contribution from the HTL photon countcrterm 
diagram with a single photon self-energy insertion is 



•^2c = ^ IP"(P)A^(P) 



(21) 



The contribution from the HTL electron counterterm 
diagram with a single electron self-energy insertion is 



2d 



{P} 



Tr [E(P)S(P)] 



(22) 



The role of the counterterm diagrams (2c) and (2d) is 
to avoid overcounting of diagrams when using effective 
propagators in (la) and (lb). Similarly, the role of coun- 
terterm diagram (3k) is to avoid overcounting when using 
effective vertices in (2a). 

The thermodynamic potential at next-to-next-to- 
lcading order (NNLO) in HTL perturbation theory can 
be written as 



^NNLO — ^NLO + N"j{T-i c + jF 3 j) + Nj{T^, a + jF 3b + Tzd + + Ftf + J-^g + + JF 3k + ^31) + JF 31l 




FIG. 4: Diagrams contributing to NNLO in HTLpt which contribute beyond order e 5 . The diagrams in the first line above 
first contribute at order e 8 and the second line at order e 6 . The undulating lines are photon propagators and the solid lines are 
fermion propagators. All propagators and vertices shown are HTL-resummed propagators and vertices. 
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2 \(dm 2 D 



2 \2 



2 1 

n LO ) {A im 2 D ) 



2 \{dm 2 ) 2 



n LO (Aim 2 ) 2 + 



F 



2a+2b 

a 



AiQ 



(23) 



where ^Sq, A2m 2 D , A2TO 2 7. and Aia are the terms of order 5 2 in the vacuum energy density, mass and coupling 
constant countcrtcrms: 
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A 2 m 2 D 



A 2 m 2 
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The contributions from the three-loop diagrams are given by 

•^3a = -^4 TY[T^-P,Q-P,Q)S(Q)r a (Q-R,Q,R)S(R)r»(P,R,R-P) 

4 JP{QB} 

xS(R - P)T f3 (-Q + R,R- P,Q- P)S(Q - P)] A^ (P)A a0 (Q — R) , 
^3b = \e 4 4 TT[T^P,P + Q,Q)S(Q)T^-R + Q,Q,R)S(R)T a (R-Q,R,Q) 

1 Tp{QR} 

xS(Q)T v (-P, Q, P + Q)S{P + Q)] A^{P)A a(3 (R - Q) , 
•^3c = -- A e*4 Tr [r^(P, P + Q, Q)S(Q)rP(-P, Q,P + Q)S(P + Q)] 

4 JP{QR} 

xTr [T u (-P, R,P + R)S(P + R)T a (P, P + R, R)S(R)} A^(P)A af3 (P) , 
Fzi = -\e iT £ Tr[T a0 (P,-P,R,R)S(R)]A a ^(P)A^(P) 

2 Tp{QR} 

xTr [r^(P, P + Q, Q)S(Q)T l/ (-P, Q,P + Q)S(P + Q)\ . 



(24) 
(25) 
(26) 



(27) 
(28) 
(29) 
(30) 



The contributions from the two-loop diagrams with -p 3f _ _t e 2 ^r Ti[T^(P —P Q Q)S(Q)] 



electron-photon three- and four-point vertices with an 
insertion of a photon self-energy 



P{Q} 

xA^ a {P)Tl af} {P)A f}v {P) , 



(32) 



•^3d = -\e 2 4 Tr[T^P,Q,R)S(Q)T"(P,Q,R)S(R)} 
1 JP{Q) 



:A^ a (P)n^(P)A^(P) 



where R = Q- P. 



The contributions from the two-loop diagrams with the 
(31) electron-photon three and four-point vertices with an in- 
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scrtion of an electron self-energy are 

^3e = -e 2 i A a P(P)Tr{T a (P,Q,R)S(Q) 

Tp(Q\ 



ip{Q} 

x^(Q)S(Q)r( j (P,Q,R)S(R)] , 



(33) 



\e 2 4 A^{P) 

2 JP{Q} 



cTr [r^(P, -P, Q, Q)S(Q)E(Q)S(Q)] ,(34) 



where R = Q-P. 

The contribution from the HTL photon countcrterm 
diagram with two photon self-energy insertions is 



■Mi 



Uf" / (p)A ua (P)u a/3 (P)A f i (P) . (35) 



The contribution from HTL electron counterterm with 
two electron self-energy insertions is 



^31 = 




Tr[S(P)S(P)5(P)E(P)] . (36) 



The remaining three-loop diagrams involving HTL cor- 
rected vertex terms are given by 



Tr 

P{Q} 



f»(P, Q, R)S(Q)T»(P, Q, R)S{R) 
xA^(P) , (37) 



^31 = -\e 2 4 Tr[T^(P,-P,Q ) Q)S(Q)} 
xA^(P) , 



(38) 



where is the HTL correction term given in Eq. (|A43[) . 
Note also that diagram (31) is the same as (2b) since there 
is no tree-level electron-photon four-vertex. 

In the remainder of the paper, we work in Landau 
gauge (£ = 0), but we emphasize that the HTL pertur- 
bation theory method of reorganization is gauge-fixing 
independent to all orders in 5 (loop expansion) by 
construction. 



scalar sum-integrals exactly was intractable and the sum- 
integrals were calculated approximately by expanding 
them in powers of mry/T and uif/T. We will follow the 
same strategy in this paper and carry out the expansion 
to high enough order to include all terms through order 
e 5 if mo and m/ are taken to be of order e. The NLO 
approximation will be perturbatively accurate to order 
e 3 and the NNLO approximation accurate to order e 5 . 

The free energy can be divided into contributions from 
hard and soft momenta. In the one-loop diagrams, the 
contributions are either hard (h) or soft (s), while at the 
two-loop level, there are hard-hard (hh) and hard-soft 
(hs) contributions. There are no soft-soft (ss) contri- 
butions since one of the loop momenta is fermionic and 
always hard. At three loops there are hard-hard-hard 
(hhh), hard- hard-soft (hhs), and hard-soft-soft (hss) 
contributions. There are no soft-soft-soft (sss) contri- 
butions, again due to the hard fermionic lines. 

In the process of the calculation we will see that there 
are many cancellations between the lower-order HTL- 
improved diagrams and the higher-order HTL-improved 
counterterm diagrams. This is by construction and is 
part of the systematic way in which HTLpt converges 
to the known perturbativc expansion. For example, one 
can see that diagrams (2c) and (3h) subtract out the 
modification of the hard gluon propagator due to the 
HTL-improvement of the propagator in diagram (la). 
Likewise, one expects cancellations to occur between di- 
agrams (lb), (2d) and (3i); (2a), (3d), (3e) and (3k); 
and (2b), (3f), (3g), and (31). Below we will explicitly 
demonstrate how these cancellations occur. 



A. One-loop sum-integrals 



Hard contribution 



IV. EXPANSION IN THE MASS PARAMETERS 

In the papers [2(| Hlj , the free energy was reduced to 
scalar sum-integrals. It was clear that evaluating these 



For hard momenta, the self-energies are suppressed by 
mo/T and rrif /T relative to the inverse free propagators, 
so we can expand in powers of Ht(P), n^P), and £(P). 

For the one-loop graph (la), we need to expand to 
second order in m 2 n : 



r l 
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-{d-\)% log(P 2 
^ jp 
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24 
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2"° D ? P P2 

C'(-i)' 
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4ttT 



(39) 



The one-loop graph with a photon self-energy insertion (2c) has an explicit factor of m 2 D and so we need only to 
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expand the sum- integral to first order in m 2 D : 



& = ^ + ^ymi J - 2^ - 2^T P + 2— T P + ^ ( T P ) 



1 

24 



2 + 2 



a-i)j \ Uttt) 



2< 



m 2 D T 2 



64tt 2 



1 n n 2tt 2 \ / u 



2c 



(40) 



The one-loop graph with two photon self-energy insertions (3h) must be expanded to zeroth order in m 2 D 



T (h) _ 

3h — 



f 



The sum of Eqs. 



4(d-l) 
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' D 



pi p 2p2 p 4 p 2p2 p 4 



f287T 2 

(|4ip is very simple 



1 „ „ 2^ 2 
- - 7 + 2 7 + — 
e 3 



4ttT 



T (h) 

- / la+2c+3h 



\{d-l)^\og(P 2 ) 



(41) 



_ _ rp4 

45 



(42) 



This is the free energy of an ideal gas of photons. 

The one-loop graph (16) needs to expanded to second order in m 2 : 
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T{p} 
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V47tT 
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12tt- 



(tt 2 - 6) 



(43) 



The one-loop fermion loop with a fcrmion self-energy insertion (2d) must be expanded to first order in m 2 , 



T. 
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Am] 
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~6 



7{ P } t> T{p } 



1+ 2-21og2 + 2 



p>4 p 2 P 2 p' 2 P 2 

C'(-i) 

C(-l) y "J ^ttT 



r P - 4^ (r P ) 2 
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(44) 



The one-loop fermion loop with two self-energy insertions 
(3i) must be expanded to zeroth order in m 2 : 
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The sum of Eqs. (f43 |) -(|45 |) is particularly simple: 
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180" 
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This is the free energy of an ideal gas of a single massless 
fermion. 



2. Soft contribution 



The soft contributions in the diagrams (la), (2c), and 
(3/i) arise from the Pq = term in the sum-integral. At 
soft momentum P = (0,p), the HTL self-energy func- 
tions reduce to Ht(P) = and Ul(P) = m 2 D . The trans- 
verse term vanishes in dimensional regularization because 
there is no momentum scale in the integral over p. Thus 
the soft contributions come from the longitudinal term 
only and read 
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(49) 



Note that we have kept the terms through order e in 
Eqs. (|47p and l|48p as they are required in the calculation 
of the counterterms. There is no soft contribution from 
the leading-order fermion term (fl"5|) or from the HTL 
counterterms (|2"2"j) and (f3"6) . 



B. Two-loop sum-integrals 



For hard momenta, the self-energies are suppressed by 
vtid/T and rrif/T relative to the inverse free propagators, 
so we can expand in powers of Ilr, I±£, and E. 



(hh) contribution 



The (hh) contribution from (fl^) and (f^Uj) was calcu- 
lated in Ref. [21J| and reads 
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(50) 

Consider next the (hh) contribution from (j3Tj) and (|32|) . The easiest way to calculate this term, is to expand the 
two- loop diagrams (2a) and (2b) to first order in ni 2 D . This yields 
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(51) 



We also need the (hh) contributions from the diagrams (3e), (3g), (3k), and (31) The first two diagrams are given 
by (|3"3"1) . (|34[) . while the last remaining ones are given by ([57)1 and (|3"8")) . The easiest way to calculate these contributions 
is to expand the two-loop diagrams (2a) and (2b) to first order in nil. This yields 
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The sum of the terms in (IST)1) -(|52 [) is very simple 
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(52) 



(53) 



This is the two-loop contribution from the perturbative 
expansion of the free energy in QED. 



2. (hs) contribution 

In the (hs) region, the momentum P is soft. The mo- 
menta Q and R are always hard. The function that multi- 
plies the soft propagator At(0, p), A^(0, p), or Ax(0, p) 
can be expanded in powers of the soft momentum p. The 
soft propagators At(0,p), Ai(0,p), and Ax(0,p) are 



defined in Eqs. ([MI]) , (TA"22f and (TA"27]) , respectively. In 
the case of At(0,p), the resulting integrals over p have 
no scale and they vanish in dimensional rcgularization. 
The integration measure J scales like mf), the soft prop- 
agators Ai(0,p) and Ax(0,p) scale like l/mjj, and ev- 
ery power of p in the numerator scales like tud- 



The terms that contribute through order e 2 m 3 D T and 
e 2 m 2 f m D T from (19]) and (20]) were calculated in Ref. HU 
and read 
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(54) 



The (hs) contribution from (|3~T]) and (32]) can again be calculated from the diagrams (2a) and (2b) by Taylor 
expanding their contribution to first order in m 2 D . This yields 
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(55) 



We also need the (hs) contributions from the diagrams (3e), (3g), (3k), and (31) Again we calculate their contribu- 
tions by expanding the two- loop diagrams (2a) and (2b) to first order in m 2 f . This yields 



q-(hs) . 2 2rp 

'3e+3g+3fc+3i ~ 4m / e 1 



pP Z + m D T{Q} 
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Q 4 Q e Q 4 Q Q 2 \(Q-Y) 2 



= ^2 m f m D T - 



(56) 



10 



3. (ss) contribution 

There are no contributions from the (ss) sector since fermionic momenta are always hard. 

C. Three-loop sum-integrals 

1. (hhh) contribution 

If all three loop momenta are hard, we can expand the propagators in powers of n^^P) and £(P). Through order 
e 5 , we can use bare propagators and vertices. The diagrams (3a), (3b), and (3c) were calculated in Refs. [2 01 and 
their contribution is 



^3aS + 3c - 2 (d 5)ei f {pQR} P 2 Q 2 R 2 (P + Q + R) 2 
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Using the expression for the sum-integrals in the Appendix, we obtain 

,57r 2 /a\ 2 • / // \ 6e 



(57) 
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(58) 



The diagrams (3a) and (3b) are both infrared finite in the limit mo — > 0. This implies that the e 5 contribution 
is given by using a dressed longitudinal propagator and bare vertices. The ring diagram (3c) is infrared divergent in 
that limit. The contribution through e 5 is obtained by expanding in powers of self-energies and vertices. Finally, the 
diagram (3j) also gives a contribution of order e 5 . Since the electron-photon four-vertex is already of order e 2 m 2 , we 
can use a dressed longitudinal propagator and bare fermion propagators as well as bare electron-photon three- vertices. 
Note that both (3c) and (3j) are proportional to N 2 and so it is more convenient to calculate their sum. One finds 
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Using the expressions for the integrals and sum-integrals listed in the Appendix, we obtain 

UttZV 
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3. (hss) contribution 

The (hss) modes first start to contribute at order e 6 , and therefore at our truncation order the (hss) contributions 
vanish. 

4. (sss) contribution 

There are no contributions from the (sss) sector since fermionic momenta are always hard. 

V. THE THERMODYNAMIC POTENTIAL 

In this section we present the final renormalizcd thermodynamic potential explicitly through order 5 2 , aka NNLO. 
The final NNLO expression is completely analytic; however, there are some numerically determined constants which 
remain in the final expressions at NLO. 

A. Leading order 

The complete expression for the leading order thermodynamic potential is given by the sum of Eqs. (|39|) . (|43|). 
and (|47[) plus the leading vacuum energy counterterm (|14p : 

^lo = |l + |iV>- f m% - 30N f m} + Wrhl + | (log \ - \ + 7 + m% - 607V, (. 2 - 6) m 4 } . 

(63) 

where mo, rhf, and /i are dimcnsionlcss variables: 

rhf = — , (65) 

B. Next-to-leading order 

The renormalization contributions at first order in 5 are 

Aifi = A 1 £ + A im l-^n LO + Aimj-^-j^o . (67) 

Using the results listed in Eqs. (|16[) . (JTTJ) , and (fT5)) . the complete contribution from the counterterm at first order in 
S is 



. _ tt 2 T 4 45 . 4 A a 

45 4e J 7r 



. (68) 



i + 2 log I - 21ogm D + 2) mi, + y Q + 2 + 2 log | - 2 log2 + 2^-^ 
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Adding the NLO counterterms (|68|) to the contributions from the various NLO diagrams, we obtain the renormalized 
NLO thermodynamic potential 
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C. Next-to-next-to- leading order 



The rcnormalization contributions at second order in <5 are 
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Using the results listed in Eqs. ([24]) . (|25| . and ([26]) . the complete contribution from the counterterms at second order 
in S is 
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Adding the NNLO counterterms (|7T|) to the contributions from the various NNLO diagrams, we obtain the renor- 
malized NNLO thermodynamic potential. We note that at NNLO all numerically determined subleading coefficients 
in e drop out and we are left with a final result which is completely analytic. The resulting NNLO thermodynamic 
potential is 
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(72) 



We note that the coupling constant counterterm listed in Eq. ([5]) coincides with the known one-loop running of 
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the QED coupling constant 
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Below we will present results as a function of e evaluated 
at the renormalization scale 2nT. Note that when the 
free energy is evaluated at a scale different than /i = 2nT 
we use Eq. ([73"]) to determine the value of the coupling at 
H = 2nT. 

We have already seen that there are several cancella- 
tions that take place algebraically, irrespective of the val- 
ues oi mo and m/. For example the (hh) contribution 
from the two-loop diagrams (2a) and (2b) cancel against 
the (hh) contribution from the diagrams (3d), (3e), (3f), 
and (3g). As long as only hard momenta are involved, 
these cancellations will always take place once the rele- 
vant sum- integrals are expanded in powers of mo/T and 
rrif /T. This is no longer the case when soft momenta are 
involved. However, further cancellations do take place if 
one chooses the weak-coupling values for the mass pa- 
rameters. For example, if one uses the weak-coupling 
value for the Debye mass, 



m 



D 




2q2 



NNLO 



(74) 



cancel alge- 



the terms proportional to 

braically and HTLpt reduces to the weak-coupling result 
for the free energy through e 5 . 



VI. FREE ENERGY 

The mass parameters mo and rrif in hard-thcrmal- 
loop perturbation theory are in principle completely arbi- 
trary. To complete a calculation, it is necessary to specify 
mo and m/ as functions of e and T. In this section we 
will consider two possible mass prescriptions in order to 
see how much the results vary given the two different as- 
sumptions. First we will consider the variational solution 
resulting from the thermodynamic potential, Eqs. © and 
(|10p . and second we will consider using the e 5 pertur- 
bative expansion of the Debye mass UH^] and the e 3 
perturbative expansion of the fermion mass (30j . 



A. Variational Debye mass 

The NLO and NNLO variational Debye mass is deter- 
mined by solving Eqs. and (flTj)) using the NLO and 
NNLO expressions for the thermodynamic potential, re- 
spectively. The free energy is then obtained by evaluating 
the NLO and NNLO thermodynamic potentials, and 
([72"]). at the solution to the gap equations ^ and (TiTj|) . 
Note that at NNLO the gap equation for the fermion 
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FIG. 5: NLO HTLpt predictions for the free energy of QED 
with Nf — 1 and the variational Debye mass. Different curves 
correspond to varying the renormalization scale \x by a factor 
of 2 around (X = 2ttT. 



mass is trivial and gives mf = 0. In Figs. [5l [6] and [7J we 
plot the NLO and NNLO HTLpt predictions for the free 
energy of QED with Nf = 1. As can be seen in Fig. [7] the 
renormalization scale variation of the results decreases 
as one goes from NLO to NNLO. This is in contrast to 
weak-coupling expansions for which the scale variation 
can increase as the truncation order is increased. 

One troublesome issue with the variational Debye mass 
is that at NNLO this prescription gives solutions for mo 
that have a small imaginary part. We plot the imag- 
inary part of the free energy which results from these 
imaginary contributions to the variational Debye mass 
in Fig. [5] (bottom panel). The imaginary contributions 
to the variational Debye mass come with both a positive 
and negative sign corresponding to the two possible solu- 
tions to the quadratic variational gap equation. The pos- 
itive sign would indicate an unstable solution while the 
negative sign would indicate a damped solution. These 
imaginary parts are most likely an artifact of the dual 
truncation at order e 5 ; however, without extending the 
truncation to higher order, it is difficult to say. They do 
not occur at NLO in HTLpt in cither QED or QCD. Wc 
note that a similar effect has also been observed in NNLO 
screened perturbation theory in scalar theories [l3[ . Be- 
cause of this complication, in the next subsection we will 
discuss a different mass prescription in order to assess 
the impact of these small imaginary parts. 



B. Perturbative Debye and fermion masses 



The perturbative Debye and fermion masses for 
have been calculated through order e 5 @, and e 3 
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FIG. 6: Real (top panel) and imaginary (bottom panel) parts 
of NNLO HTLpt predictions for the free energy of QED with 
Nf = 1 and the variational Debye mass. Different curves 
correspond to varying the renormalization scale by a factor 
of 2 around jj, = 2nT. 



FIG. 7: A comparison of the renormalization scale variations 
between NLO and NNLO HTLpt predictions for the free en- 
ergy of QED with Nf = 1 and the variational Debye mass. 
The bands correspond to varying the renormalization scale fi 
by a factor of 2 around fi = 2nT. 
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Plugging dxSj> and (J76J) into the NLO and NNLO thermo- 
dynamic potentials, ([55)) and {72]), we obtain the results 
shown in Fig. [5J The renormalization scale variation is 
quite small in the NNLO result. 



FIG. 8: A comparison of the renormalization scale variations 
between NLO and NNLO HTLpt predictions for the free en- 
ergy of QED with Nf = 1 using the perturbative thermal 
masses given in Eqs. (|75p and (|76p . The bands correspond to 
varying the renormalization scale p by a factor of 2 around 
H = 2ttT. 



C. Comparison with the ^-derivable approach 

Having obtained the NNLO HTLpt result for the free 
energy we can now compare the results obtained using 
this reorganization with results obtained within the 4>- 
derivable approach. In Fig. [9] we show a comparison 
of our NNLO HTLpt results with a three-loop calcu- 
lation obtained previously using a truncated three-loop 
<&-derivable approximation [3l|. For the NNLO HTLpt 
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FIG. 9: A comparison of the predictions for the free energy 
of QED with Nf = 1 between three-loop "^-derivable approx- 
imation [H and NNLO HTLpt at [i = 2nT. 



FIG. 10: A comparison of the predictions for the free energy 
of QED with Nf = 1 between the two- loop 2PI approximation 
in Landau gauge [H and NLO HTLpt at n = 2vrT. 



prediction we show the results obtained using both the 
variational and perturbative mass prescriptions. As can 
be seen from this figure, there is very good agreement 
between the NNLO ^-derivable and HTLpt approaches 
out to large coupling. In all cases we have chosen the 
renormalization scale to be = 2ttT. 

As a further consistency check, in Fig. [TU] we show a 
comparison between the untruncated two-loop numeri- 
cal ^-derivable approach calculation of Ref. [32j and our 
NLO HTLpt result using the variational mass. In both 
cases we have chosen the renormalization scale to be 
/i = 2ttT. From this figure we see that there is a reason- 
able agreement between the NLO numerical ^-derivable 
and NLO HTLpt results; however, the agreement is not 
as good as the corresponding NNLO results. 

We note that the results of [32| were computed in the 
Landau gauge (£ = 0). As detailed in their paper, their 
result is gauge dependent. Such gauge dependence is un- 
avoidable in the 2PI ^-derivable approach since it only 
uses dressed propagators. In Ref. [3l| it was explicitly 
shown that the two-loop ^-derivable Debye mass is gauge 
independent only up to order e 2 , resulting in gauge vari- 
ation of the free energy at order e 4 . This is in agreement 
with general theorems stating that the gauge variance 
appears at one order higher than the truncation (33j . 



VII. CONCLUSIONS 

In this paper we calculated the three-loop HTLpt ther- 
modynamic potential in QED. Having obtained this we 
applied two mass prescriptions, variational and perturba- 
tive, to fix the a priori undetermined parameters mo and 
rrif that appear in the HTL-improved Lagrangian. We 
found that the resulting expressions for the free energy 



were the same to an accuracy of 0.6% at e = 2.4 giving 
us confidence in the prediction. We also compared the 
HTLpt three-loop result with a three-loop ^-derivable 
approach [3l| and found agreement at the subpercentage 
level at large coupling. 

In addition, we showed that the HTLpt NLO and 
NNLO approximations have improved convergence at 
large coupling compared to the naively truncated weak- 
coupling expansion and that the renormalization scale 
variation at NNLO using both the variational and per- 
turbative mass prescriptions was quite small. Therefore, 
the NNLO HTLpt method result seems to be quite re- 
liable. This is important since, unlike the ^-derivable 
approach, the HTLpt reorganization is gauge invariant 
by construction and is formulated directly in Minkowski 
space allowing it to, in principle, also be applied to the 
calculation of dynamical quantities. 

The renormalization of the three-loop thermodynamic 
potential required only known vacuum, mass, and cou- 
pling constant counterterms, and the resulting running 
coupling was found to coincide with the QED one-loop 
running. This provides further evidence that the HTLpt 
framework is rcnormalizablc despite the new divergences 
which are introduced during HTL improvement. 

Finally, we note that at three loops we could obtain an 
entirely analytic expression for the renormalized NNLO 
thermodynamic potential. There were a number of can- 
cellations that took place during renormalization which 
resulted in an expression that was independent of any nu- 
merically determined subleading coefficients in the sum- 
integrals. We expect similar cancellations to also occur 
in non-Abclian gauge theories which will greatly simplify 
the calculation. Computing the three-loop HTLpt reor- 
ganized free energy for QCD is in progress. 
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APPENDIX A: HTL FEYNMAN RULES 

In this appendix, we present Fcynman rules for HTL 
perturbation theory in QED. We give explicit expres- 
sions for the propagators and for the electron-photon 
three- and four-vertices. The Feynman rules are given 
in Minkowski space to facilitate applications to real-time 
processes. A Minkowski momentum is denoted p = 
(po, p), and the inner product is p ■ q = poqg — p • q. The 
vector that specifies the thermal rest frame is n = (1, 0). 



1. Photon self-energy 

The HTL photon self-energy tensor for a photon of 
momentum p is 



(Al) 



The tensor T fiu (p, q), which is defined only for momenta 
that satisfy p + q = 0, is 



T^(p,-p) 



p-y 



(A2) 



The angular brackets indicate averaging over the spatial 
directions of the lightlike vector y = (l,y). The tensor 
is symmetric in /i and v and satisfies the "Ward 
identity" 



(A3) 



The self-energy tensor Ii^ v is therefore also symmetric in 
fi and v and satisfies 

(A4) 
(A5) 



Pa JP» = 0, 
g^Wip) = -rr 



The photon self-energy tensor can be expressed in 
terms of two scalar functions, the transverse and lon- 
gitudinal self-energies and Hl, defined by 

u T ( P ) = - fp>) ir j ( P ) , (A6) 

H L (p) = -H 00 (p) , (A7) 



where p is the unit vector in the direction of p. In terms 
of these functions, the self-energy tensor is 

IP"(p) = -IItW - \n L (j>)L^ , (A8) 



where the tensors T p and L p are 



P 



p^p 



iP"n v 



The four-vector np 1 is 



(A9) 
(A10) 

(All) 



and satisfies p-n p = and = 1 — (n-p) 2 /p 2 . Equa- 
tion (|A5|) reduces to the identity 



1 



(d-i)n r (p) + -2n L (p) = m 



(A12) 



We can express both self-energy functions in terms of the 
function T 00 defined by |X2|) : 

™° [T m (p,-p)-l + nl] , (A13) 

(A14) 



(d- l)n| 
U L ( P ) = m 2 D [l-T 00 (p,-p)}, 



In the tensor T M!/ (p, —p) defined in (|A2[) , the angular 
brackets indicate the angular average over the unit vec- 
tor y. In almost all previous work, the angular average 
in (|A2|) has been taken in d = 3 dimensions. For consis- 
tency of higher-order radiative corrections, it is essential 
to take the angular average in d = 3 — 2e dimensions and 
analytically continue to d = 3 only after all poles in e 
have been cancelled. Expressing the angular average as 
an integral over the cosine of an angle, the expression for 
the 00 component of the tensor is 



2 J-i Po 
where the weight function w(e) is 



PC 



r 2 (i 



rmr(i-e) 



(A15) 



(A16) 



The integral in (|A15|) must be defined so that it is ana- 
lytic at po = 00 • It then has a branch cut running from 
po = — |p| to po = +|p|- If we take the limit e — > 0, it 
reduces to 



r°°(p,-p) = ^iog^ 

2|p| Po — |p| 



(A17) 



which is the expression that appears in the usual HTL 
self-energy functions. 
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The Feynman rule for the photon propagator is 

iA^Cp) , (A18) 

where the photon propagator tensor A^„ depends on the 
choice of gauge fixing. We consider two possibilities that 
introduce an arbitrary gauge parameter £ : general covari- 
ant gauge and general Coulomb gauge. In both cases, the 
inverse propagator reduces in the limit £ — > oo to 



p^p y -W u {p). (A19) 



This can also be written 



/ 1 V 



A t (p) 



n 2 p A L (p) P 



K v > ( A2 °) 



where At and Al are the transverse and longitudinal 
propagators: 



A T Op) 
Al(p) 



1 



p 2 - Tl T (p) 
1 



-n 2 p 2 



n L (p) ' 

The inverse propagator for general £ is 
A" 1 ^)^ = A" 1 ^)^ - ^p»p v covariant, 



(A21) 
(A22) 

(A23) 



p-n n") 
Coulomb . (A24) 



The propagators obtained by inverting the tensors 
in (|A"54l) and (|A"23|) are 



A""(p) 



-At{p)T£ u + Ai(p)n^rip — £ 



p^p" 



(P 2 ) 2 
covariant , (A25) 



= -A T (p)T j r + A i (p)n'V'-e 



p>y 
n 2 p 2 ) 



Coulomb . (A26) 

It is convenient to define the following combination of 
propagators: 

A x (p) = A L {p) + \a t { p ) . (A27) 



Using (|A12|) . (|A21|) . and (|A22|) . it can be expressed in 
the alternative form 

A x (p) = [m 2 D - dU T {p)] Ai(p)A T (p) , (A28) 

which shows that it vanishes in the limit mo — > 0. In the 
covariant gauge, the propagator tensor can be written 

A""(p) = [-A T (p)g^ + A x (p)nW] 

-^A x (p)(pV + #/) 



P 



At(p) + ^#A x (p)-1 
p z p^ 



(A29) 



This decomposition of the propagator into three terms 
has proved to be particularly convenient for explicit cal- 
culations. For example, the first term satisfies the iden- 
tity 

[-A T (p)g^ + A x (p)n^}A x 1 ( P r x = 
a P^P A . n-p A x (p) A 



3m 



p 2 n 2 p 2 A L (p) 



Pfl n; . (A30) 



2. Electron self-energy 

The HTL self-energy of an electron with momentum p 
is given by 



S(P) = m 2 r(p) , 



where 



T"(p) 



P ' 2/ 



(A31) 



(A32) 



Expressing the angular average as an integral over the 
cosine of an angle, the expression is 



T"(p) 



w(e) f 1 



Po - |p|c 



(A33) 



The integral in (|A33[) must be defined so that it is ana- 
lytic at po = oo. It then has a branch cut running from 
po = — |p| topo = +|p|- In three dimensions, this reduces 
to 



2 P Po - p 



/ „ Po , Po + P 

-r47 ■ P 1 - ^i-r log H 

P V 2P po- P 



(A34) 



3. Electron propagator 

The Feynman rule for the electron propagator is 

iS{p) . (A35) 
The electron propagator can be written as 



Sip) 



1 



i - s(p) 



(A36) 



where the electron self-energy is given by (|A31|) . The 
inverse electron propagator can be written as 



5" 1 (p)=^-S(p) 
This can be written as 

S~ 1 (p)=4(p), 



(A37) 



(A38) 
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where we have organized Aq(p) and As(p) into: 

A li (p) = (A (p),A s (p)i>) . (A39) 
The functions Ao(p) and A sip) arc defined 



m* 

A (p) = Po L T P , 

Po 



Mp) = | P | + ^[i-r p ] 



(A40) 
(A41) 



4. Electron-photon vertex 

The electron-photon vertex with outgoing photon mo- 
mentum p, incoming electron momentum q, outgoing 
electron momentum r, and Lorentz index n is 

r"(p, q, r) = e (> - m 2 f»{ Pl q, r)) . (A42) 

The tensor in the HTL correction term is only defined 
for p — q + r = 0: 



t"(p > «,r) = (^ 



(A43) 



This tensor is even under the permutation of q and r. It 
satisfies the "Ward identity" 



p^( P ,q,r)=T»(q)-T»(r) 



(A44) 



The electron-photon vertex therefore satisfies the Ward 
identity 

p^(p,q,r) = S- 1 (q)-S- 1 (r) . (A45) 



5. Electron-photon four- vertex 

We define the electron-photon four-point vertex with 
outgoing photon momenta p and q, incoming electron 
momentum r, and outgoing electron momentum s. It 
reads 

T^(p,q,r,s) = -2e 2 m 2 f^(p 7 q,r,s) 

= 2e 2 r^, (A46) 

There is no tree-level term. The tensor in the HTL cor- 
rection term is only defined for p + q — r + s = 0, 



Y[ — + — 

v-y s-y 



[{r-p)-y] [(s+p)-y] 



.(A47) 



This tensor is symmetric in (i and v and is traceless. It 
satisfies the Ward identity: 

p^ v (p, q, r, a) =T v (q, r-p,s)- V v {q, r, s + p) .(A48) 



6. HTL electron counterterm 

The Feynman rule for the insertion of an HTL electron 
counterterm into an electron propagator is 



i£(p) 



(A49) 



where is the HTL electron self-energy given 

in |[A39"1) . 

7. Imaginary-time formalism 

In the imaginary-time formalism. Minkowski energies 
have discrete imaginary values po = i(2TrnT) and inte- 
grals over Minkowski space are replaced by sum-integrals 
over Euclidean vectors (27rnT, p). We will use the nota- 
tion P = (PoiP) for Euclidean momenta. The magni- 
tude of the spatial momentum will be denoted p = |p|, 
and should not be confused with a Minkowski vector. 
The inner product of two Euclidean vectors is P ■ Q = 
PqQo + P " Q- The vector that specifies the thermal rest 
frame remains n = (1,0). 

The Feynman rules for Minkowski space given above 
can be easily adapted to Euclidean space. The Euclidean 
tensor in a given Feynman rule is obtained from the corre- 
sponding Minkowski tensor with raised indices by replac- 
ing each Minkowski energy po by iPo, where Po is the cor- 
responding Euclidean energy, and multipying by —i for 
every index. This prescription transforms p = (po,p) 
into P = (Po,p), into -5^, and p-q into -P-Q . 
The effect on the HTL tensors defined in (jA2|, i[A43]) . 
and (|A47[) is equivalent to substituting p-n —> —P-N 
where N = (— i,0), p-y — > —P-Y where Y = (— i, y), 
and y^ — > Y^. For example, the Euclidean tensor corre- 
sponding to (|A2j) is 



/ P-N 
T^iP, -P) = (YW'-py 



(A50) 



The average is taken over the directions of the unit vector 

y- 

Alternatively, one can calculate a diagram by using 
the Feynman rules for Minkowski momenta, reducing the 
expressions for diagrams to scalars, and then make the 
appropriate substitutions, such as p 2 — > — P 2 , p ■ q — > 
— P • Q, and n ■ p — > in ■ P. For example, the propagator 
functions (|A21|) and (|A22|) become 



A t(P) = ^ 



-1 



Al(P) = 



P 2 + n T (p) ' 
i 

P 2 + n L (P) ' 



(A51) 
(A52) 



The expressions for the HTL self-energy functions Hr(P) 
and U L (P) are given by (|A13|) and (|A14|) with n 2 re- 
placed by n 2 p — p 2 /P 2 and T°°(p, —p) replaced by 



W{€) f l dc (1 - c 2 )- lP ° 



iP - pc 



(A53) 
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Note that this function differs by a sign from the 00 com- 
ponent of the Euclidean tensor corresponding to (|A2|) : 



T UU (P,-P) = -T uu (p, -p) 



-T P 



(A54) 



po— nPo 



A more convenient form for calculating sum-integrals 
that involve the function Tp is 



T P = 



P 2 



P Q 2 + p 2 c 2 



(A55) 



where the angular brackets represent an average over c 
defined by 



4 ; 

T{P } 



eV 

47T 



P = (2n+l)irT 



2< , 



(27T) 



3~2£ 



(B2) 



where 3 — 2e is the dimension of space and is an arbi- 
trary momentum scale. The factor (e 7 /47r) € is introduced 
so that, after minimal subtraction of the poles in e due 
to ultraviolet divergences, p, coincides with the renormal- 
ization scale of the MS renormalization scheme. 

Below we list the sum-integrals required to complete 
the three-loop calculation. We refer to [10, HH for details 
concerning the sum-integral evaluations. 



(f(c)) c = w(e) / dc(l-c 2 r/(c), 
Jo 

and w(e) is given in (|A1GI) . 



(A56) 



1. One-loop sum-integrals 



APPENDIX B: SUM-INTEGRALS 

In the imaginary-time formalism for thermal field the- 
ory, the four- momentum P = (Po,p) is Euclidean with 



P 2 
o 



p . The Euclidean energy po has discrete 



p2 

values: Po = 2ri7rT for bosons and Po = + 1)ttP for 
fermions, where n is an integer. Loop diagrams involve 
sums over Po and integrals over p. With dimensional 
regularization, the integral is generalized to d = 3 — 2e 
spatial dimensions. We define the dimensionally regular- 
ized sum-integral by 



47T 



T 



E 

Pa=2mrT 



d 3 ~ 



P 



(27T) 



3-2£ 



(Bl) 



The simple one-loop sum-integrals required in our cal- 
culations can be derived from the formulas 



H \2e 2T(| +m- e)r(n - § - m + e) 



(P 2 )™ 



4ttTV r(n)r(2-2e) 
x T(l - e)C(2n - 2m - 3 + 2e)e £7 
x T 4+2m - 2 "(27r) 1+2m - 2 " , (B3) 



< { P}(P 2 ) n (2 Tp (P 2 ) n ' 

The specific bosonic one-loop sum-integrals needed are 



(B4) 



logP 2 = 



p 

fpP 2 
1 

pW? 
i 



-^T 4 [l + 0( e )] 



T / fi 



12 V4ttT 



2+2 cHy) 



(4tt) 2 V47rT 



2c 



- + 2 1 + 0{e) 



1 



/' 



, P p 2 P 2 (An) 2 \4irT 
The specific fermionic one-loop sum-integrals needed are 

7n 2 



- + 2 7 + 2 + C(e; 
e 



t 

T{p} 



\ogP z 



1 

P}P~ 2 
1 



360 

T 2 



T 4 [1 + 0(e)] 



/' 



24 \47tT 



1+ 2-21og2 + 2 



C'(-i) 
C(-i) 



{P} 



(47T) 



— f— V 



- + 2 7 + 41og2 + C(e) 
e 



0( 



(B5) 

(B6) 

(B7) 
(B8) 

(B9) 
(BIO) 
(BH) 
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P~ 



{P} (P 2 ) 2 
P 2 



11 
16 

1 



f— ) 



2< 



1 



2^3 



{P} (P 2 ) 
P 4 



(4tt) 2 4 
/" \ 2c 5 



1 



'{p> 



i 



(4tt) 2 V4ttT/ 8 
1 



p 2 P 2 (4tt) 2 V4ttT 



4 C'(-l) s 



- + 2 7 - § +41og2 + C(e) 

- + 2 7 --^+41og2 + 0(e) 
e 15 

- + 2 + 2 7 + 41og2 + C(e) 



0(e 



(B12) 
(B13) 
(B14) 
(B15) 



The errors are all of one order higher in e than the smallest term shown. The number 7l is the first Stieltjes gamma 
constant defined by the equation 



1 



C(l + «) = -+7-7i* + O0O 



(B16) 



We also need some more difficult one- loop sum-integrals that involve the HTL function defined in (|A33[) . The 
specific bosonic sum-integrals needed are 



p P 



p 2p2 
1 



T P = 



1 fJJ_ 
(4tt) 2 UttT 

(4tt) 2 UttT 



-1) 



- + 2 7 + 21og2 + C(e) 



1 



1 



/' 



r P p- 1 ' (4tt) 2 UttIV V 3 

The specific fcrmionic sum-integrals needed are 



( ^ + 2 7 j +21og 2 2 + y +<D(e) 

(1 \ 4 22 
(l + 21og2) (- + 2 7 J - - + ylog2 + 21og 2 2 + 0(e) 



1 



{P} (P 2 ) 2 

f 1 

{P}^ 

i T 



(4tt) 2 UttT/ 2 

1 / A* N2e 

(4tt) 2 UttT 

1 / /i 



(4tt) 2 V4ttT 



{ p } P 2 ^o 2 

X P} P 2 \(p-r) 2 / 



(4tt) 2 V4ttT 
(4tt) 2 UttT 



- + 2 7 +l + 41og2 + C(e) 
e 



21og2 ( - + 2 7 J +101og 2 2+y +C( e ; 



1 1 7T 2 

+ 2( 7 + 21og2)- + — +41og 2 2 + 8 7 log2-4 7 i + 0(e) 
e z e 4 

2 , 



log2 [ - + 2 7 ] +51og z 2 + C(e) 



) 



- - l + 2 7 + 41og2 + C(e) 
e 



(B17) 
(B18) 
(B19) 

(B20) 
(B21) 

, (B22) 
(B23) 
(B24) 



2. Two- loop sum-integrals 

The simple two-loop sum-integrals that are needed 
r i 



are 



= o, 



?{pq} ^ 2 Q 2 r 2 (4tt) 2 WtJ 

f g 2 t 2 ( » \ Ae ( M 

{PQ} ^ 2 Q 2 ^ 4 ' ' (4tt) 2 UttT ) \ 12 J 



£ T 2 / M \ 4e / 1 

{PQ} P 2 Q 2 r 2 R 2 (4^) 2 UttT/ ^ 72 
P-Q T 2 



f JL 



{PQ} P 2 Q 2 r 4 (4^) 2 UttT 



4f 



■ +4 -21og2 + i^—j + 0(e) 
1 11 



y + 2 7 -21og2 + 2^^+0(e) 
7.00164 + 0(e) 



I-67 + 6^^ + 0(6) 



(B25) 
(B26) 
(B27) 
(B28) 
(B29) 
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T 



{PQ} q 2 P 2 Q 2 R 2 (4tt) 2 \AttTJ \72 



{PQ} Q 2 P 2 Q 2 R 2 (4tt) 



— m 

'4tt) 2 V47tT/ 



4( 



1 
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9.55216 + C(e) 



8.14234 + 0(e) 



(B30) 
(B31) 



where R = — (P + Q) and r = |p + q|. The corrections are all of order e 2 . We also need some more difficult two-loop 
sum- integrals that involve the functions Tp defined in (|A33[) . 



1 



{ p Q} 



P 2Q2 r 2 
,2 



q 



{PQ} P'Q'r 4 

PQ 

P 2Q2 r i 



{PQ} 



2 2 



{PQ } PWQIR 2 



Tr 



Tr 



Tr 



tt) 2 UttTV 



4r 



(4 

T 2 
(4tt) 2 

T 2 



1) 


1 


487 


e 2 " 


— ) 


1 


576/ 


e 2 


— ) 


1 


96/ 


e 2 ~ 


1 





2 + 12 log 2 + 4 -^-jy ) ^ + 136.362 + 0(e) 
^ + 52 log 2 + ) - + 446.412 + 0(e) 



41og2 + 4^-^ ) i +69.174 + 0(e) 

C(-i)y e 



10 



C(-i)A 1 



2 + 2 7+y log2 + 2 c( _ i) 



(4tt) 2 UttT/ 

3. Three-loop sum-integrals 



+ 46.8757+ 0(e) 



(B32) 
,(B33) 

(B34) 
. (B35) 



The three-loop sum-integrals needed are 



i 



1 


P 2 Q 2 i? 2 (iM 
1 


Q- 


hi?) 2 


P 2 Q 2 R 2 (P^ 


Q- 


hi?) 2 


1 






P 2 Q 2 R 2 (P^ 


Q- 


hi?) 2 


QR 







{p }Q R P 2 Q 2 R 2 {P + Q) 2 {P + Rf 



(Q-Rf 



p{qr } P 2 Q 2 R 2 {P + Q) 2 {P + Rf 



(47 



-f-Vr-^v 

tt) 2 I 12 } UttT/ 



1 / T 



(4tt) 2 V 12 7 V4ttT 

,C'(-i) 



,(B36) 



+12 



C(-i) 



0(e) 



(B37) 



1 f T W 



(4tt) 2 V 12 7 V4ttT 
+0(e 

i f T2 Y ( M 



3 179 51 3 C(-3) C C'(-1) 

i e -lb- + io log2+ 2cF3T~ 6 cRy 



(4tt) 2 V 12 7 V4ttT 
3C'(-1) 



3 9 361 57, „ 3C'(-3) 
^ + 4^ + M) + T0 1Og2+ 2CF3T 



2 C(-l) 



0(e) 



1 



/ M \ 6£ 



(4tt) 2 V 12 y V47rT7 
, 7C'(-1) 



5 

24e + 4 



1 23 8 1CC-3) 
7 + — — — log 2 — 



24 5 



6 C(-3) 



6 C(-l) 



0(e) 



(B38) 



(B39) 



(B40) 



The three-loop sum-integrals were first calculated by Arnold and Zhai, and calculational details can be found in 
Rcf. 0. 



APPENDIX C: THREE-DIMENSIONAL INTEGRALS 



Dimensional regularization can be used to regularize both the ultraviolet divergences and infrared divergences in 
three-dimensional integrals over momenta. The spatial dimension is generalized to d = 3 — 2e dimensions. Integrals are 
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evaluated at a value of d, for which they converge, and then analytically continued to d = 3. We use the integration 
measure 



eV 

47T 



'P 



(27T) 



3-2t 



(CI) 



The one-loop integrals needed are of the form 

In. 



Specifically, we need 



P (p 2 + m 2 ) 



2\n 



1 

8^ 



(e^V) e 



T(n - | + e) 

r(|)r(n) 



,3-2n-2e 



(C2) 



I' = [ log(p 2 + m 2 



h 
h 



TO 3 / fl ^ 2e 



6ir \2m 
m / /i 

1 / fX x2 



l + -e + 0(e 2 ) 



■f f [l + 2 £ + 0(e 2 )] 
47T \2m/ L v /J 



87rm 



(/-) [l + 0(e)] 
\2mJ 
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